Dielectric response due to stochastic motion of pinned domain walls 
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We study the contribution of stochastic motion of a domain wall (DW) to the dielectric AC 
susceptibility for low frequencies. Using the concept of waiting time distributions, which is related 
to the energy landscape of the DW in a disordered medium, we derive the power-law behavior of 
the complex susceptibility observed recently in some ferroelectrics below Curie temperature. 
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During the last decade considerable interest has 
been attracted to the anomalously high dielectric 
susceptibility observed in several ferroelectrics below 
the paraelectric-ferroelectric phase transition temper- 
ature Tc- The theory of critical phenomena pre- 
dicts in the vicinity of Tc the Curie- Weiss law for 
the temperature dependent static susceptibility. How- 
ever, for some ferroelectrics, such as potassium di- 
hydrogen phosphatsi (KH2PO4) or rubidium dihydro- 
gen phosphate^ (RbH2P04), the complex susceptibil- 
ity x{^) = x'('^) ~ *x"(^) within the so-called plateau 
range Tf < T < Tc remains unusually high. This 
anomaly is usually ascribed to DW motion, an ac- 
tivity which freezes out below Tf. According to 
Ref. Ill, the low-frequency dielectric spectrum in the 
plateau range of KH2PO4 consists of non-Debye contri- 
bution, which has a power-law behavior with a small 
exponent, and several Debye-like constituents. Re- 
cent experiments performed on RbFei/2Nbi/203'^ and 
Sro.6i-xCexBao.39Nb206 {SBN:Ge)— revealed in the 
range of low frequency lo < lJc a. power-law behavior of 
the DW response 



XH=Xoo(l + l/(«C.Tor) 



(1) 



with < 71 < 1. Here lOc is the threshold frequency 
separating the low- and high- frequency regimes, tq the 
characteristic relaxation time and Xoo the high-frequency 
limit of the complex susceptibility. The real and imag- 
inary parts of the complex susceptibility (^ are related 
by equation 



X"M = (x'H - Xoo) tan(nV2), 



(2) 



which follows from the Kramers-Kronig relation. In a 
range of frequencies above u>c, the imaginary part of 
susceptibility x"{'^) ^^^ found to increase linearly on a 
log-log scale, while the corresponding real part x'('^) de- 
creases linearly on a linear-log scale 1^ 

An attempt to explain phenomenologically the exis- 
tence of two regimes of the dynamic response was made 
in Ref. ,5, where the low-frequency behavior (^ was at- 
tributed to the irreversible viscous motion of DWs in a 
disordered medium. However, this interpretation, based 
on the picture of the dynamic hysteresis developed in 
Ref. la, is restricted to the adiabatic regime only. 



Eqs. (E3) and Q are in accordance with Jonscher's uni- 
versal dielectric response lawji^ which describes the re- 
sponse of a wide class of dielectrics in a broad range of fre- 
quencies. Generally, this behavior is related to hopping 
of charge carriers, such as electrons, polarons or various 
ions, between different localized states. The exponent n 
is determined by the spatial and energetic distribution 
of the corresponding states. The ferroelectric non-Debye 
dielectric response ^ is expected to be due to jumps 
of DWs (or DW segments) between different metastable 
states. 

In the present paper, we study the contribution of 
stochastic motion of a DW to the dielectric susceptibil- 
ity at low frequencies, using the concept of waiting time 
distribution, which is related to the energy landscape of 
the DW in a disordered medium. We restrict our con- 
sideration to the regime in which the field driven motion 
of adjacent DWs is independent of each other. Accord- 
ing to Ref. y the exponent n depends on the poling state 
of dielectric samples. Poling is expected to decrease the 
density of DWs, which would improve the single DW ap- 
proach. However, among others it may also change the 
distribution of impurities. 

The DW, which separates domains of different polar- 
ization, can be viewed as a two-dimensional (d — 2) elas- 
tic interface moving in a disordered environment. The 
configuration of a d-dimensional elastic interface is de- 
scribed by the profile function z{x,t), which obeys the 
Langevin-type equation 

^-^ij;^^2_)_ = j^^z + hoCOSLut -\- g{x, z) -^ 'q{x,t), (3) 

where fi is the mobility, 7 the stiffness constant and 
h{t) = ho cos cut the AC driving force density. Note that 
X is the d-dimensional vector. We assume Gaussian dis- 
tribution of thermal noise, with zero mean and the cor- 
relator 



r]{x,t)r]{x',t') = 2fi-^TS^'^\x ~ x')S{t - t'). (4) 

Note that in real ferroelectrics the dependence of /i and 
7 on temperature T may complicate the comparison of 
the computed temperature dependence of the suscepti- 
bility with experiment. Pinning to impurities, which are 
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FIG. 1: The energy landscape of the center of mass of DW 

described by the quenched random force g{x,z), strongly 
suppresses the DW mobility. The quenched forces are 
assumed to be Gaussian distributed with zero mean and 
the variance 



{g(x, z)g{x', z')) = S'-'^^ {x - x')A{z - z'). 



(5) 



To make this model well-defined, we introduce the cutoff 
Aq in the Dirac's delta function 5'^{x) at scales in the 
order of the impurity separation. In what follows we 
consider random field (RF) disorder, where for z > the 
correlator A(z) — A(— z) is a monotonically decreasing 
function of z, and decays rapidly to zero over the distance 
a. 

To describe different regimes of the DW motion, we 
recall the fundamental scales of the problem. On scales 
smaller than the Larkin scale Lc — (a^7^/A(0))^/^, the 
elasticity wins over the disorder for d < 4 and the DW re- 
mains flat, while on larger scales the DW becomes rough 
and self-affine with a nontrivial roughness exponent C,.^'^ 
At T = and ui ^ 0, the system ©-© undergoes a 
second order depinning transition at ho = /i^iiiiiSii^ The 
depinning threshold he — ajL~'^ can be estimated from 
balancing the pinning forces and the driving forces on the 
scale Lc- At finite temperature, the depinning transition 
smears out, and the DW moves at any field strength. In 
the creep regime, which sets in for low fields, the motion 
is very slow and controlled by thermal activation. The 
energy barrier Eb{L) ~ Uc{L/Lc)^ has to be overcome 
in order to move a DW segment of linear size L, where 
9 = 2( + d - 2, and Uc ^ 0^7^^"^ is the typical bar- 
rier on the Larkin scale Lc- The DW segment of size L 
spends time t{L) = tq exp(£'s(i)/T) in a valley before 
it jumps to the next valley. Here tq = Ll/{iJ,j) is the mi- 
croscopic time. For low fields h <$: he, the average DW 
velocity is determined by the thermal activation of seg- 
ments with optimal size iopt = ic(^/^c)~^ , and 
obeys the creep law 



vo,eM-Ue/T{h/her% 



(6) 



with/i= (d- 2 + 20/(2 -C). 

If the DW is driven by AC field with the frequency 
uj, the barriers with waiting times u)t{L) > 1 cannot 
be overcome during one period of the field. Thus the 
threshold frequency 

to* = l/r(iopt) = r^^ eM~Uc/T{ho/he)-n (7) 



separates the sliding regime from the pinned regimei^ 
The critical value uj*, which corresponds to the exper- 
imental setup used in Refs. and I3, can be estimated 
as follows. Since the activity of DWs for T < Tf is 
almost frozen, ^^ we conclude that temperatures in the 
range Tf < T < Tc are of order Uc- The creep ex- 
ponent is iJ,{d = 2) = 1 for RF disorder. The ampli- 
tude of the electric field in the dielectric experiments 
of Ref. ^ho~ 300V/m is far below the coercive field 
150kV/m. Assuming that the threshold he is of order 
of the coercive field we find that the experimental fre- 
quencies (w > O.OlHz) are likely to be much higher than 
the threshold frequency estimated from Eq.(|7I). To ex- 
plain the observed low-frequency response we suppose 
that DWs can also move for lu* < uj < cUe {^^c is defined 
below) via hopping between different minima of the en- 
ergy landscape. However, this motion is extremely slow 
and can be characterized by the displacement z oc i" with 
n < 1^ so that in the adiabatic limit t ^ 00 the aver- 
age velocity will vanish. In analogy with the stochastic 
transport phenomena in disordered solids j— we refer to 
the DW-hopping in the frequency range to* < uj < LUe SiS 
the stochastic regime. 

To calculate the dielectric response in the stochastic 
regime, we consider a 180° domain structure of quasiperi- 
odicity 21, where the DW separates homogeneously polar- 
ized regions with spontaneous polarization Pq and —Po, 
respectively. The time dependence of the macroscopic 
polarization reads 



dP{t)/dt^{2Po/l){v{t)). 



(8) 



The linear susceptibility x{^) defined through the Fourier 
transform P{uj) — x{^)h{uj) of P{t) can be expressed as 



x{^) 



2Po fijuj) 
I iuj 



(9) 



where fi{uj) is the renormalized DW mobility given by 
{v{uj)) — fi{uj)h{uj)- The center of mass of the field driven 
DW probes different local minima of the rugged energy 
landscape, corresponding to different metastable DW- 
configurations (see Fig. QJ . This motion is similar to the 
stochastic motion of a particle in the RF environment^ 
The Langevin equation for the particle is given by 



dz 
'di 



f(z) + h + rjit), 



(10) 



where f{z) is the Gaussian distributed RF with {f{z)) — 
0, {f{z)f{z')) — AqS{z — z'), and ri{t) the thermal noise. 
At low temperatures T and small driving forces h, the 
motion of the particle can be described by the thermally 
activated dynamics, which is controlled by the probabil- 
ity distribution density W{E) of energy barriers E sep- 
arating different metastable states. Using the standard 
methods of the field theory, W{E) can be written as a 
path integral, the estimate of which by the steepest de- 
scent method gives^^ 



W{E) 



Eq ^ exp( 



-E/Eo) 



(11) 



with Eq = Ao/{2h). Due to the thermally activated dy- 
namics, the waiting time t{E) = Toexp{—E/T) needed 
to overcome the barrier E is characterized by the distri- 
bution 'i>{T)dT = W{E)dE, which decays as a power law 
for large r 



*(r) = 



T-l+n' 



n - T/Eo 



(12) 



As follows from Eq. H12|l the average waiting time 
J cItt'^It) diverges for low temperatures T < Eq. This 
indicates that the particle motion is dominated by very 
rare, but extremely deep potential wells, which leads to 
sublinear drift z ex t", n < 1. 

The distribution of waiting times which controls the 
motion of the center of mass of DW can also be related to 
the distribution of the energy barriers W{E)ii^ As shown 
in Ref. 18, the energy barriers scale as the energy minima, 
and therefore, have identical distributions. Generalizing 
the loop expansion of the energy distribution developed 
in Ref. 19 for driven elastic interface we have obtained the 
energy distribution of the DW in the presence of nonzero 
external force h as 



WiE) 



ds 

2^ 



exp(zsi?) exp 



-II^ 



X I In 

Ik 



IS 

7fc2 



1-f g2/j2/(^2fc4) 



,(13) 



where Ag is the Fourier transform of A(z). Assum- 
ing that the width a of the disorder correlator Aq is 
small, we approximate the integral over q in Eq. H13|) 
by A(0)/(1 -t- h?a^ /^"^k^). The inverse Fourier trans- 
form of the discrete version of Eq. (|13(l for d = 1 
gives f{x]K) = YlT=i ^A>^)'^^V{-{P + i^'^/P)x), where 
X = Att'^-iE/{/\{Q)L'^), k ^ {haL'^)/{'iTr^-f), and the coef- 
ficients Cj{k) are given by 



C,{k) 



2j 



2j-5 



r(2j)r J - - r 



X if - .-) ^°'^ (" 



,c2 



)-cos2(7ryf) 



1=1 



(14) 



\f 



Using the steepest descent method to estimate the sum 
in f(x;K), we find that for large energy E the distribu- 
tion is exponential W{E) oc exp(— 2a/ii?/A(0)). How- 
ever, this expression of the energy distribution is incor- 
rect for the same reason as the perturbation result for 
the roughness?12. The appropriate method to extract the 
correct behavior of elastic interfaces in random environ- 
ment is the functional renormalization group iS^iiiSiiiiS 
To find the effect of renormalization on the energy distri- 
bution for driven interface, we follow Ref. [l^ and replace 
all bare quantities in Eq. 1)131) by the renormalized ones. 
We obtain the renormalized energy distribution for d — 1 
as W{E) oc exp(-2(7a/i)i/2£'/A(0)Lc)- The evaluation 
of Eq. (|13|l for d > 1 is more complicated. However, for 



d = 2 the tail of the energy distribution for large energies 
is exponential, and is given by Eq. (|ll|l with Eq of order 
Uc- Consequently, the distribution of waiting times has 
the same form (|12f) . 

To calculate the renormalized mobility ^{lo), we 
now employ the MontroU- Weiss formalism developed in 
Refs. 16 and \2'A for continuous time random walks. Ac- 
cording to Ref. Ila. the mobility /x(a;) can be written as 






dtt 



'^{{z{t)-z{Q)f) 



2Tl-§(zu;)' 



(15) 



where uq is the average distance between adjacent poten- 
tial wells and \['(u) is the Laplace transform of ^(t). For 
^(r) given by Eq. H12|l we obtain 



^(u) = / dte""^*(T) ~ 1 - A(uTo)", 



(16) 



where A — 7r/(r(ri) sin(7rn)). Combining Eqs. 0, 115)1 . 
and l)16)l . we obtain the complex susceptibility as 



x(w) 



2Poal 



Al 2T {iuT^y- 



(17) 



The above consideration applies for low frequencies, 
where the center of mass of the DW explores many val- 
leys during one period 2tt /to . For uj larger than some 
frequency uJc, the center of mass of the DW stays in one 
valley. The characteristic frequency tUc can be estimated 
from the condition that the center of mass of the DW re- 
mains within one well during period 2tt/lu with the prob- 
ability 1. This yields lUc ^ Tq^ . 

For frequencies oj > lUc, the DW is captured in one 
valley, the potential of which we approximate by 



"^'>'w/ 



dz'g{z'), 



(18) 



where the 1st term describes the average shape of the 
well, and the 2nd term describes the fluctuations due to 
the residual part of the random force g. Neglecting the 
second term in Ea. l)18)) yields Debye-response 



X"(w) oc UJTl/{l + (uJTlf). 



(19) 



The 2nd term in Eq. ()18l) results in distribution of re- 
laxation times ^i(t) instead of the single relaxation time 
Ti. Expecting the exponentially distributed barriers on 
all scales, the distribution of relaxation times is likely to 
obey a power-law ^i(t) — niT"^/r^+"i with some expo- 
nent Hi which is not necessarily equal to n. As a result 
the Debye response ()19)l changes to 



X"(a;) a (c.ri)"^ 



(20) 
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FIG. 2: The imaginary part of the susceptibility of poled 
SBN:Ce measured at T = 294 Ki^ The best fit by Eq. ^ 
corresponds to a = 21.56, n = 0.69, k = (n + n\)/2n = 0.75 
(ni = 0.35), and ujc ~ 91. 



To describe x"(w) for w ^ Wc we propose the following 
scaling ansatz 



x"{io) = {LUT,r^F{Ljao). 



(21) 



F{x) is a scaling function and L^ is the distance in z 
direction explored by the DW during one period of the 
AC field. The distance L{Tya) explored by the DW during 
the time Tm can be estimated from the condition 



L{t^) 



cro 



'^{T)dT = 1, 



(22) 



which means that the DW at most once is captured by 
the deepest well with the waiting time larger than t^. 
For Tm equal to the period of AC field, Ea. lf^ gives 
i(j = i(Tni ~ l/oj) = ffQ / (lutq)" . Substituting the latter 
into Eq. 121|l we obtain 



x"{iu) = {cjT,r^F{i/{cjTor). 



(23) 



According to Eas. H17|l and H2U|) . the scaling function F{x) 
must behave for small and large x as 



Fix) 



x"S a;<l, 
j,i+«i/n^ a;> 1. 



The fit of the experimental data of Kleemann et al.^ 
using the scaling function 



x"{o.) 



aoj 



1 + 



2k- 1 



(24) 



is shown in Fig. El 

In contrast to the imaginary part, the real part of the 
susceptibility behaves as a power of Inw.-^ To explain 



this behavior, we suggest that the intrawell dynamics of 
DW is a superposition of center of mass motion and relax- 
ational motion of internal modes. Indeed, in the pinned 
phase uj > lOc, the DW segments can also jump between 
metastable states with close energies, which gives rise to 
additional dissipation. The relaxational contribution of 
internal modes to the susceptibility reada^S^ 



X' ex (ln(l/c.ro))'/' , 
X" « (ln(l/c.ro))'/'-^ . 



(25) 
(26) 



We briefly remind the derivation of Eqs. H25|l and H26|l . 
As suggested in Ref s . |2l| and ^(f. one can treat the pinned 
DW as an hierarchical ensemble of noninteracting two- 
level systems (TLS), each is a DW segment of linear 
size L. The separation between the two configurations 
of the given TLS is w{L) ex L^, and the energy dif- 
ference is AE. Transitions caused by thermal activa- 
tion may occur only in TLSs with the energy difference 
AE < T. The corresponding rate of transitions is given 
by T-\L) = Tf^^exp{-EB{L)/T), where Eb{L) is the 
energy barrier. In the presence of AC field, these transi- 
tions lead to power dissipationSi 



Q{L) 



{5E{L)f 



lo^t{L) 



4Tcosh2(A/2r)l+wV2(L)' 



(27) 



where SE = PQL'^w{L)ho is the shift of AE due to the ex- 
ternal field h. To obtain the power dissipation density, we 
have to average Eq. (|27|) over the distribution P£(A_E), 
and then sum over all TSLs in the unit volume. Since the 
distribution function Pl{AE) is expected to be smooth 
with a width of the order AE{L) (x L^ ^ T, the mean 
value of the factor cosh~^(A£^/2T) gives the fraction of 
the thermally active TLSs T/AE{L). The hierarchical 
structure of TLSs implies that the density of TLSs of 
linear size L is 1/L''+^. Therefore the power dissipation 
density Q can be written as 



Q oc Pohl / dLL 



LU\{L) 

'1 + cjV2(L)' 



(28) 



Due to the exponential increase of r with L, the main 
contribution to the integral (|28|l comes from the length 
scale L(j given by the condition ujt(Li^) ^ 1. On the 
other hand, the power dissipation is related to the imag- 
inary part of the complex susceptibility^^ 



7=f " J. 2 



(29) 



Equating (|28|) and H29|l we arrive at Eq. (|26|l for the imag- 
inary part of the susceptibility. The corresponding real 
part H25I) can be restored from the imaginary part l|2(j|) 
with the help of the Kramers-Kronig relation in the form 
of"7r/2rule"i2i 



X"{u;) 



2 din a; 



(30) 



The real part of the susceptibility, which corresponds 
to Eq. (|20|l . decreases as a power-law and can be hidden 
by the slower decreasing logarithmic contribution H25|l . 
In contrast, the power-law contribution 12U|) dominates 
the behavior of the imaginary part of the susceptibility, 
so that the logarithmic contribution (|26(l is irrelevant. 

Let us finally discuss in brief the nonlinear response 
associated with the DW motion. We define the nonlinear 
susceptibilities as Xm = x'm ~ *Xm with 



1 



2Tr 



X'm = ^j P{t) cos{mujt)d{ujt), (31) 



^™ = ^ / ^(^) Mrniot)diLot). (32) 



The linear susceptibility corresponds to xi- To calculate 
the nonlinear susceptibility in the sliding regime, i.e. at 
very low frequencies uj < uj*, we integrate Eq. (jSj) with 
(v(i)) given by Eq. © and h{t) = hQCOs{ujt), and then 
substitute P{t) in Eqs. l|^ and lj221).^^' As a result we 
obtain x'm — ^'^d Xm °^ ^/^ for odd m. Note that 



in the linear case (m = 1) this gives a conduction- type 
susceptibility xi'^) — ^/iujTQ. The computation of the 
nonlinear susceptibility in the stochastic regime is a more 
difficult task which is left for further investigations. 

In conclusion, we have considered the dielectric re- 
sponse of ferroelectric DWs below the paraelectric- 
ferroelectric phase transition temperature. As a function 
of the frequency of the external electric AC-field there 
are three regimes in the behaviour of the complex sus- 
ceptibility: the sliding, the stochastic, and the pinned 
regimes. The response in the sliding regime, which cor- 
responds to very low frequencies, is due to the creep- like 
motion of DWs. The response in the stochastic regime is 
related to jumps of the DW as a whole, and can be de- 
scribed using the concept of waiting time distributions. 
The response in the pinned regime at high-frequencies 
is expected to be due to superposition of center of mass 
motion and relaxational motion of internal modes. The 
existence of all three regimes for the dielectric response 
was recently observed in SrTi Os^ 
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